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Abstract 



Q \ Motivated by the effect of matter on neutrino oscihations (the MSW effect) 

we study in more detail the propagation of neutrinos in a dense medium. 
^ ' The dispersion relation for massive neutrinos in a medium is known to have a 

minimum at nonzero momentum p ~ Cpp/y^. We study in detail the origin 
and consequences of this dispersion relation for both Dirac and Majorana 
neutrinos both in a toy model with only neutral currents and a single neutrino 
flavour and in a realistic "Standard Model" with two neutrino flavours. We 
find that for a range of neutrino momenta near the minimum of the dispersion 
relation, Dirac neutrinos are trapped by their coherent interactions with the 
medium. This effect does not lead to the trapping of Majorana neutrinos. 
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I. INTRODUCTION 



Motivated by the effect of matter on neutrino oscillations (the MSW effect |T]]), there 
have been several works in recent years aimed at understanding in a more complete way 
the propagation of one or more flavours of massive neutrinos in matter. One of the first 
papers along these lines was the paper of Mannheim in 1987 whose main purpose was 
to derive the MSW effect from a Field Theoretic starting point. Mannheim used second 
quantization techniques to derive the wave functions and dispersion relations of two flavours 
of both Dirac and Majorana neutrinos propagating in a medium in which there was a finite 
density of electrons. Mannheim then analyzed his result in the ultrarelativistic regime and 
recovered the standard MSW results. 

Both the work of Mannheim and the work of Nieves and of Notzold and Raffelt 0] 
showed that the entire MSW effect could be reliably analyzed with a modified Dirac (or 
Majorana) equation by adding to the frequency of the electron neutrino a term proportional 
to the density {\/2GfP in the standard MSW scenario). This term is analogous to a chemical 
potential term for the electron neutrino though it couples only to the chiral left handed part 
of a possibly massive neutrino. 

In 1991 Panteleone [^] used such a modified Dirac equation to study the behaviour 
of neutrinos in supernova cores. He first analyzed the case of a single neutrino flavour 
without restricting the momenta to be relativistic. At this point he noticed a very unusual 
characteristic of the neutrino dispersion relation. The dispersion relation had a minimum 
at a nonzero neutrino momentum. Thus for a range of neutrino momentum the neutrino's 
phase velocity and group velocity are in opposite directions! Panteleone later analyzed the 
case of two and three flavours of neutrinos. He noticed that the neutral currents did not 
decouple in general although his detailed analysis was carried out only in the high energy 
regime in which they do decouple. 

In this paper we provide a synthesis of many of the above results. We are particularly 
interested in effects at low neutrino momentum and in effects due to the minimum of the 
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dispersion relation at nonzero momentum. We begin in Section |T| by examining a simple 
model with only a single neutrino flavour in which the neutrino propagates in a background 
of electrons. We take the neutrino-electron interaction to be mediated only by neutral 
current interactions. This model allows for a careful analysis of the Field Theory aspects 
of neutrino propagation in a medium. We pay special attention to the minimum of the 
dispersion relation which occurs at non-zero momentum and analyze its effects on neutrino 
interactions. We will find that the minimum energy for Dirac neutrinos in the medium will 
generically be less than the rest mass of the neutrino in the vacuum so that very low energy 
neutrinos are effectively trapped by the medium. A similar effect has been studied from a 
different point of view in the work of Loeb [] 0, although his effect has an entirely different 
origin. Our analysis of the trapping of neutrinos will bring up many interesting questions 
and puzzles which will need to be resolved in order to have a complete understanding of 
the problem. We will also examine the case of Majorana neutrinos and find that in general 
Majorana neutrinos are not trapped. 

In Sec. ^ we extend our analysis to a more realistic model which has two neutrino 
flavours and in which there are both neutral and charged current interactions. In this case 
the dispersion relations are governed by quartic equations for Dirac neutrinos and quadratic 
equations for Majorana neutrinos. We will again be most interested in examining the form of 
the dispersion relations in the medium and their effects on neutrino propagation. Again we 
will find that the Dirac case leads quite generically to neutrino trapping while the Majorana 
case can have no trapping. We will also make a few remarks regarding the oscillations of 
neutrinos in a medium, noting that even for a single neutrino there could in principle be 
oscillations in the probability to detect the neutrino, due to the differing phase velocities of 

^Loeb studies the problem from the point of view of the neutrino's "index of refraction" in the 
medium. A radially varying density leads to a force on the neutrino which allows it to have bound 
orbits in the medium. 
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the helicity eigenstates. 

We conclude in Sec. |V| with a brief discussion of our results and some concluding remarks. 



II. A SIMPLE MODEL WITH ONE NEUTRINO FLAVOUR 

Many of the interesting effects which we will discuss, namely those due to the minimum 
of the dispersion relation which occurs at non-zero momentum, occur already in a simple 
model with only a single neutrino flavour. It is useful, then, to first consider a simplified 
model in which a Dirac neutrino propagates in an electron "gas" to which it couples only 
via the neutral current interaction. The case of a Majorana neutrino is somewhat more 
subtle and will be considered subsequently. Our model may be described by the following 
Lagrangian 

£ = (z^ - m,) + i'e - me) ' + - /ieV'eVe, (2.1) 

where 

Dt = ^,±^^Z,il-^% (2.2) 
i^M- = d^Z, - d,Z^. (2.3) 

The chemical potential term in the Lagrangian is included in order to give a non-zero value 
to the electron density; that is 

Pe = (^eVe) 7^ 0. (2.4) 

In order to study the propagation of a neutrino in this medium, we compute the neutrino 
self-energy. To one loop there are three diagrams, shown in Fig. |l| 0. All effects due to the 

^The electron loop in this diagram is, of course, equal to the electron density to all orders. The 
neutrino loop and the Zq loop lead to a term proportional to the neutrino density which, as we 
have pointed out, is nonzero. This leads to a small correction which can easily be included but 
which we shall ignore. 
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non-zero electron density come from the electron loop in Fig. |l|(a) which is easily calculated 
and yields 



S = -^Pe7°(l-7'), 



(2.5) 



in which we have defined G=g^ j (4-\/2w,|) in analogy with the usual Fermi coupling constant, 
Gp. From the self-energy one may obtain the neutrino propagator in the usual way by 
summing a geometric series. For constant pe the resulting expression is given in momentum 
space by 

1 



(2.6) 



J - - S 

If pe depends explicitly on x the propagator may still be formally written in position space 
as 

1 



GAx) 



i(f) — m — S(a;) 
The effective action is then given, to this order, by 



eff 



G-Hx] 



G 



^+^Pe7°(l-7') 



(2.7) 

(2.8) 
(2.9) 



Variation of the effective action leads finally to an effective "Dirac equation," given by 

"z^-m + a7°(l -7^)] ^/-^ = 0, (2.10) 

in which we have defined a=Gpe/ V2. For constant electron density, the presence of the 
"chiral potential" in this expression leads to a shift in the frequency by a, but only for the 
left-handed (chiral) piece. This shift in the frequency is precisely the "index of refraction" 
familiar from the MSW effect. Once we have derived the dispersion relations, it will be 
clear that the shift in energy for the neutrino is opposite that for the anti-neutrino 0. If the 
neutrino is "attracted" by the medium, then the anti-neutrino is "repelled" by it. 



^This shift comes from a term in the effective Hamiltonian proportional to iplipu which equals the 
number density of neutrinos minus the number density of anti-neutrinos. 
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We have noted above that the chemical potential /le in the Lagrangian (|2.1| ) gives rise 
to a non-zero electron density, Pf,=((ip\%ljf,) . A similar calculation of py={iplilj^) for the ef- 
fective Lagrangian defined by Eq. ( p.9[ ) shows that, at least naively, there appears also to 
be a non-zero density of neutrinos in this medium (provided, of course, that the effective 
chemical potential a is larger than m). It is clear that this arises due to the potential term 
proportional to 7*^ in the effective Lagrangian, since this term looks exactly like a chemi- 
cal potential for the chiral left-handed neutrinos. How one handles this apparent density 
of neutrinos can drastically affect the MSW effect. Suppose, for example, that we were to 
insist that in the sun Pu={'ipl'ipiy)=0. In order to implement this, we may choose to introduce 
a "counter" chemical potential into the original Lagrangian which would exactly cancel the 
neutrino density generated by the interactions with the electrons in the medium. It turns 
out, however, that this would change the MSW result quite dramatically. In fact, if our 
theory had a "vector" instead of a "chiral" potential, we would kill the entire MSW effect 
by doing this. As we shall discuss in more detail below, the correct thing to do is to accept 
the fact that in equilibrium there is a non-zero density of neutrinos of very low momentum 
due to their attractive interaction with the medium [^j. 

Further insight into the physics of our model can be gained by examining the equations 
of motion following from the Lagrangian in Eq. (|2.1| ). Varying the Lagrangian with respect 
to the Z'^ field leads to 



d^F^f" + mlZ'' = ^Jf", (2.11) 

where 

= ^,^{1 - 7')^e - V^.7^(l - 7')^.- (2.12) 
If Pe={4'l4'e) is constant, then ( |2.11| ) leads to 



that is, the Z^ field has gained a vacuum expectation value. The equation of motion for the 
neutrino field is then given by 



iq) — m — 



2^/2 



(ZV(l-7') 



= 0, 



(2.14) 



which is equivalent to the effective Dirac equation of ( ^.1U| ) once the value for the expec- 
tation value, Eq. ( |2.13| ), is inserted. From this point of view, then, the left-handed neutrino 
sees a mean (coherent) "scalar potential," (Z^). From the point of view of the field theoretic 
calculation above, it is clear why the Z^ field has developed a vacuum expectation value. 
The one- loop diagram corresponding to (Z^) is simply the electron loop, corresponding to pe, 
with a Z propagator attached. This coherent Z° field is analogous to the electric field which 
surrounds a static charge distribution and is due to the net weak charge of the medium. 

A. Solution to the Dirac Equation 

In order to study the propagation of neutrinos over macroscopic distances, it suffices to 
study the effective Dirac equation given in Eq. ( |2.10| ). The propagator defined in Eq. ( f^.7| ) 
contains this information, but it also encodes the off-shell behaviour of the neutrino. For a 
medium with constant density, it is straightforward to solve the Dirac equation in momentum 
space by employing the chiral representation, so that 



[xr) 



(2.15) 



in which the upper and lower components correspond to the left and right chiral projections, 
respectively. In this representation, the Dirac equation becomes 



— m uj — a ■ p 
uj + 2a + a ■ p —m 



0. 



(2.16) 



\XrJ 



where a=Gpe/ V2. Without loss of generality we may choose p=pz, so that Xl,r (and hence 
also tp) may be chosen to be eigenstates of cx^, the spin projection in the z direction. That 



IS, 



Cr3XL,R = SXl,R, 



(2.17) 



where s=±l. Solving for the energy yields four solutions 

uj = -a±^{p + asy + m^. (2.18) 

These dispersion relations are plotted in Fig. ^ both for m=0 (dashed curves) and rriy^O 
(solid curves.) Several key features of these plots should be noted. First of all, the "negative 
energy" states are, in this case, those which are unbounded from below as the momentum 
is increased. In the second quantized theory the correct energy of such a state is just the 
negative of its energy eigenvalue. We also note that when m=0 there are "level crossings." 
These are avoided for m^O by level repulsion due to the mixing of the levels. 

The most noteworthy feature of these dispersion relations (discussed previously by Pan- 
taleone) is the fact that the minima of the dispersion relations occur at non-zero values of 
the momentum, p=±a, instead of at the origin. One interesting consequence of this fact is 
that the neutrino can have a vanishing group velocity at non-zero momenta. Furthermore, 
for |p|<|q;|, the neutrino's group velocity, duj/dp, is in a direction opposite to its momentum! 
This will play an important role in understanding the reflection of neutrinos at the bound- 
ary of the medium. Another interesting feature of these curves is that the minimum energy 
|^min=~cfc + m is Icss than the neutrino mass. Thus it is possible to produce a neutrino in 
the medium which has uj<m. Such a neutrino will not have enough energy to survive in the 
vacuum and will thus be trapped by the medium. We shall examine these peculiar features 
of the neutrino dispersion relations in detail below. 

Finally, we note that for high momentum, the solution corresponding to negative helicity 
has energy 

2 

uj^p+ 2a, (2.19) 

2p 

which is just the usual MSW result. By way of contrast, the positive helicity solution 
approaches its vacuum value of ujK.p + rn? /2p. This illustrates the spin-dependence of the 
interaction. For high momentum, the left-handed (chiral) states are nearly equivalent to the 
negative helicity eigenstates, so the potential (which is left-handed) affects only the negative, 
and not the positive, helicity eigenstates. 
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B. Neutrino Trapping 



Let us now consider a low energy neutrino which is produced inside the medium and 
which then tries to escape to the vacuum. We take the medium to have a>0. The 
dispersion relations for this neutrino both inside and outside the medium are illustrated 
in Fig. ^ In this figure the neutrino energy as a function of its momentum is plotted for 
zero density (outside the medium) (Fig. ^(a)), low density (Fig. 0(b)), intermediate density 
(Fig. |^(c)) and high density (Fig. |^(d)). The negative energy solutions to the Dirac Equation 
are also plotted, for reasons which will become apparent momentarily. 

Let us begin by looking at the case of low electron density, 0<a<m. This region is 
characterized by the fact that the minimum of the neutrino dispersion relation lies above 
the vacuum energy {E=0) but below the neutrino mass {E=m). In this case a neutrino 
which is produced inside the medium with momentum p such that E{j))>m will escape from 
the medium (although there will also be some amplitude for refiection). If, on the other 
hand, its momentum p is close to a so that E{p)<m (as is shown, for example, by the point 
A in Fig. ^b)), then this neutrino is trapped in the medium since there is no Energy level 
corresponding to its energy if it escapes into the vacuum (Fig. 0(a)). The way this total 
refiection is realized in practice is fascinating. Suppose the neutrino is incident normal to 
the interface of the medium and the vacuum. We might then be concerned that the neutrino 
has to fiip its spin in order to find an energy level with negative momentum even in the 
medium (see Fig. It is straightforward to show, however, that the neutrino's spin must 
be conserved in such a reaction. The resolution of this apparent problem is that the refiected 
neutrino indeed has a positive momentum as shown by the point B in Fig. ^(b). The reason 
this corresponds to a refiected wave is that the neutrino's group velocity is negative in this 
region so that the neutrino travels back into the medium. It is straightforward to extend 
these arguments to the case of non-normal incidence 0. 



^Note that for a<0 the medium would trap anti-neutrinos instead of neutrinos. 
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The case of intermediate electron density, m<a<2m, is shown in Fig. ^(c) and is char- 
acterized by the fact that the minimum neutrino energy in the medium is less than zero 
but greater than the largest negative energy state outside the medium. In this case the 
lowest energy of the system occurs when all the neutrino states below E=0 are filled so 
that the mean neutrino density p^, is nonzero. This is consistent with the field theoretic 
calculation of this density which was discussed previously in this paper. In this case the 
neutrinos with m>E>0 would be trapped as in the previous case but no neutrinos with 
E<0 could be produced in the medium since all corresponding levels are full. One could 
imagine a nonequilibrium (higher energy) situation in which these levels are not full. In this 
case neutrinos with energy E<0 (due to their attractive interaction with the medium) could 
be produced (by some external reaction) and they would also be trapped. Whether such 
a situation arises in any practical case depends on the dynamics of the formation of this 
region. 

Finally we turn to the case of high electron density for which a>2m. This situation 
is shown in Fig. ^(d) and is characterized by the fact that the minimum of the neutrino 
dispersion relation inside the medium has a lower energy than the maximum energy of the 
negative energy Dirac states in the vacuum. In this case the potential difference between the 
vacuum and the medium is strong enough to induce pair production of neutrino-antineutrino 
pairs at the interface. If the neutrino levels below E=—m are not filled initially then this 
pair production would eventually lead to the filling of these levels. If a neutrino with energy 
below — m is produced by another mechanism during this time, it would be trapped by the 
medium. Q 

In all the cases discussed above it is simple to write down the criterion which must be 



^In this case the reflection amplitude would also have a contribution from pair production which 
would allow a spin flip at the boundary |^. This contribution will be the dominant contribution 
in the case of massless neutrinos. 
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met in order that a neutrino be trapped. The dispersion relation for a negative spin (relative 
to z) neutrino is given by 

00= -a + ^{p-ay + m^, (2.20) 
so that the condition for trapping is 

-a + i/(p-a)2 + m2 < m (2.21) 

or 

p < ptrap = a + \Ja^ + 2am. (2.22) 

Thus, neutrinos produced in this medium with momentum p<Ptrap will not have enough 
energy to survive in the vacuum and will be trapped. 

Before proceeding it is useful to get some idea of the overall magnitude of the effect of 
neutrino trapping. Setting G^Gp and m~10~^eV (which is a mass relevant for the MSW- 
resolution of the solar neutrino problem), we find that ptrap~10^^eV in the sun (for which 
a~10^^2eV) and ptrap~100eV in a supernova (for which a~100eV.) The phenomenon of 
trapping is quite remarkable when we consider that the mean free path (which increases 
with decreasing momentum) of a neutrino with p~10^^eV in the sun is on the order of 10^° 
solar radii. Such a neutrino would thus have no chance of being "incoherently" trapped in the 
sun (say by back-scattering from nuclei), but would still be trapped by the coherent process 
which we are discussing. We note furthermore that in the case of incoherent trapping the 
scattering cross section is typically dependent on the mass of the target particle, whereas for 
coherent trapping this is not the case. Loeb has discussed a similar effect (albeit of different 
origin) for neutron stars and has estimated that in general this effect will add an extra 30 
kg to the mass of the star ||^. 

C. The Majorana Case 

In the case of Majorana neutrinos, there is only a single (left-handed) field, xl- The 
effective Lagrangian is given by 

11 



(2.23) 



where 



XL 



(2.24) 



V 



in the chiral representation. In general one needs to be somewhat careful when dealing with 
Majorana fermions. For example, even at the classical level, the fields need to be taken as 
Grassman-valued, or else the mass term disappears. The dispersion relations in this case 
can be obtained by solving the equations of motion, as was first done by Mannheim [Q. The 
reader is referred to Mannheim's paper for details of the calculation (see also Ref. |jlO|.) The 
resulting expression for the negative helicity neutrino is given by 



In this case the energy has a minimum value, 00=171, which occurs at |p|=2a. In fact, the 
dispersion relation in matter is identical to that in vacuum except for a lateral shift to the 
right. This implies in particular that, in contradistinction to the Dirac neutrino 
cannot have an energy less than m and there is thus no trapping of Majorana neutrinos in 
the medium. 

III. DISPERSION RELATIONS FOR TWO NEUTRINO FLAVOURS 

We turn now to consider a more realistic scenario in which there are two neutrino fiavours 
and in which there are both neutral current and charged current couplings to the medium. 
We have in mind, of course, the Standard Electroweak Model with massive neutrinos. We 
shall first derive the quartic equation governing the dispersion relations in the Dirac case 
and examine the solutions in some representative cases. We shall see again that in this case 
there is neutrino trapping. We then examine the Majorana case, in which the dispersion 
relations are quadratic and are thus readily analyzed but in which both trapping and pair 
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production at the surface of the sun are absent. We shall then describe briefly an alternate 
model which has been studied in the literature which has no chiral coupling but still yields 
a minimum in the dispersion relation at non-zero momentum. Finally, we comment on 
neutrino oscillations in these models. 



A. Dirac Case 



We begin with the Dirac Equation in the mass basis for a pair of massive Dirac neutrinos 
with both neutral and charged current coupling to a medium: 

{^-M+iP- aQ) (1 - 75)} V; = (3.1) 

where M is the diagonal 2x2 mass matrix, (3 oc pGp is the contribution of the neutral 

current which couples only to the left handed neutrinos and a oc PeGp represents the charged 

current contribution which couples only to Vf.- This coupling is assured by the mixing matrix 

/ cos2(0) sin(e) cos(^) \ 
Q=\ ■ (3.2) 

Vsin(0)cos(0) sin2(0) / 

a and (3 may be calculated by computing the one loop contributions to the neutrino self- 
energy in the background medium. The Feynman diagrams corresponding to these processes 
are shown in Fig. | and yield [§,||,|ll| 

a = (3.3) 

and 

/? = E(^P - 2Q^f^sm'ew)Pf (3.4) 

in which the sum in ( p.4| ) runs over all fermions in the medium, T3 is the third component 
of the fermion's weak isospin and Q^^^ is its charge. If there are appreciable densities of 
anti-particles in the medium, then pf needs to be replaced by Pf—pj in these expressions. 
In the Chiral representation we write 
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^ = ( 1 (3.5) 
\Xr. 

and the Dirac Equation becomes 

{uj-a-p)xR = MxL, {ico + a-p) + 2iP-aQ)}xL = MxR. (3.6) 

For simplicity we may assume the momentum to be in the z direction in which case the 
solutions to the Dirac Equation will be eigenstates of a^: 

leading to the equations: 

{u^ -p^ + 2{u TP){P- aQ) - M^} L± = 0, (3.8) 

R± = -^—ML±. (3.9) 



To find the energy eigenvalues we rewrite Eq. ( p.8| ) as: 

[uj^ - - fi^ + 2{uj T P)P - 2a{uj T P)N^} L± = (3.10) 

where /x^ = (m^) = {ml + 1712) /2 is the mean squared mass and 

/cos2(^)-^^ sin(e) cos(^) \ 

Nt={ , (3.11) 
Vsin(^)cos(e) sin2(^)+^^/ 

with 

= T-^ ^ (3-12) 

and = — ml. It thus remains only to find the eigenvalues of A^=p. 
The eigenvalues of A^^^ are: 

= ^ (1 + - 4^T (cos(2e) - ^^)) (3.13) 



= - ( 1 - Vl - (cos(20) - e^) ) • (3.14) 



14 



Plugging these back into the Dirac Equation leads to the following quartic equation for the 
energy eigenvalues: 

cu^ - + (2/3 - a){uj - sp) 

in which s=±l is the eigenvalue of a^, the spin projection in the +z direction. In special 
cases this expression reduces to those found in the papers of Mannheim (in which the neutral 
current contribution has been left out) and Pantaleone (in which one of the masses has been 
set to zero.) 

Equation (|3.15|) is the main quartic equation whose eight solutions (for s=±l) lead to 
the eight dispersion relations (four positive energy and four negative energy) in the medium. 
It may not at first be obvious that all eight of the solutions u of ( p.l5|) corresponding to 
a fixed value of the momentum are real. This is, however, the case, which may be seen as 
follows. (The proof is equally straightforward for any number of flavours, so we will do it 
immediately in the general case.) In the general case the Dirac Equation in the mass basis 
becomes 

- M + (/? - aU^NM) 7° (1 - 75) } ^ = (3. 16) 

in which U is the mixing matrix in flavour space and A'^e=diag(l, 0, . . . , 0) is also a matrix 
in flavour space. Pre-multiplying this expression by 7° leads to the eigenvalue equation 

Mi) = u^, (3.17) 

where 

TV" = 7°7 ■ p + -f'^M -[(3- aU^NM) (1 - 7^)- (3.18) 

Since M is hermitian for real p, the eigenvalues of Eq. ( |3.17| ) are guaranteed to be real. This 
completes the proof. 

It is best to analyze the expression governing the dispersion relations, Eq. ( |3.15| ), by first 
considering some special cases. The very simplest case is when mi=m2=0, which yields 

15 



= a^{uj - spf - aA^ cos{29){cu - sp) + -A^ (3.15) 



u = p+{p-a){s-l) 



(3.19) 



uj = -p- {p-a){s+ 1) 



(3.20) 



for the electron neutrinos and 



u = p + P{s-1) 



(3.21) 



u = -p ~ p{s + 1) 



(3.22) 



for the muon neutrinos. These expressions are easy to understand. Four of the dispersion 
relations are unchanged from their values in the vacuum, since positive helicity neutrinos 
are also right-handed (chiral) and are thus unaffected by the left-handed Standard Model 
interactions. The remaining four dispersion relations are displaced vertically from their 
vacuum values by amounts proportional to their couplings to the medium. Note that only 
the dispersion relation corresponding to Ue is affected by the charged current contribution, 
a. 

Another simple case occurs when the coupling 6 is set to zero. In this case the dispersion 
relations for z/g and z/^ decouple, as one would expect, and we find 



for the electron and muon neutrinos, respectively. These expressions are in exact agreement 
with what we found in the single-neutrino case in Sec. [II A| . Once again the dispersion 
relations corresponding to the massless case undergo "level repulsion" when a finite mass is 
added. 

Since the equations governing the dispersion relations in the two-flavour Dirac case are 
quartic it is difficult, in general, to obtain analytic expressions for these dispersion relations. 
Of course quartic equations are analytically solvable and we know that in our case all the 




(3.23) 



and 




(3.24) 
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solutions are real. In general, however, no practical insight can be gained by examining the 
analytic expressions of these solutions. One approach which is helpful in understanding the 
dispersion relations if the coupling 9 is not too large is to use a graphical approach. One 
begins by looking at the solutions when ^=0 in which case the two neutrino flavours decouple 
and we can use the solutions derived in the previous section. Thus, for example, in Fig. |^(a) 
the dotted curves represent the solutions for ^=0. Note that the dispersion relations for the 
two flavours of neutrinos cross at some points. When 9 is "turned on" we expect that these 
levels will repel and will lead to a curve similar to the solid curve in that figure. The solid 
curve is, in fact, the solution to the quartic equation when ^=0.2. This graphical method is 
reasonably accurate when 9 is small but can be used as a guide even for larger values of 9. 

One interesting feature in the two-neutrino case is that the heavier neutrino, which would 
have been trapped in certain cases if the neutrinos were decoupled, can now "leak out" due 
to its coupling to the lighter mass eigenstate. That is, only states with energy less than the 
mass of the lightest mass eigenstate are strictly "trapped" now. 

It is also possible to derive approximate solutions of the quartic equations if the neutrinos 
are relativistic. In that case approximate solutions are given by 

,2 1 r , , o -|l/2 



u;~p-(2/3-«) + ^^±-^ 



(Aap - cos(2^)) + sm^{29) 



4ap + A^ cos(2^) ) + ASin2(2^) 



1/2 



(3.25) 
(3.26) 



where the corrections to the above expressions go like a/x^/p^, and /i^/p^. Of these 

expressions, ( p.25| ) gives the energy of the negative-hehcity particle eigenstates and (p.26|) 
gives the energy of the positive-helicity antz-particle eigenstates. These are in agreement with 
the usual result and show that the neutral current contribution "factorizes" in the relativistic 
limit; that is, the difference between the two negative- helicity particle energies is independent 
of 13. Furthermore, it is clear that, if a;>0 (which occurs if the background contains more 
electrons than positrons) then a resonance can occur when 4aj9=A^ cos(20). This is the 
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well-known MSW resonance. The remaining dispersion relations, given in Eq. ( |3.27|) , are 
unchanged from their vacuum values (since the potential is left-handed) and are related to 
the positive-helicity neutrinos and negative-helicity anti-neutrinos. 



B. Major ana Case 



The case of Majorana neutrinos is interesting for two reasons. First of all, it is the 
favoured realistic scenario in models which have massive neutrinos, for example in models 
which employ the "see-saw" mechanism. Secondly, it turns out that the equations governing 
the dispersion relations are quadratic rather than quartic, which means that in principle they 
should be easier to analyze. 

The calculation proceeds in a manner similar to that followed in Sec. |II C| , the only 
complication being the additional mixing in flavour space. We omit the details and simply 
present the result. The negative-helicity dispersion relations in this case are determined by 
the equation 



where 



{u'-p'-Alip)) {u'-p'-Al{p))=0, 



AKp) = ^(ml + ml + Ap{a - 2/3) + Aa{a - 2(3) + 8/3^ 



±- 



[ml - ml) cos(26') - iap - 4a(a - 2(3) 



+ (777-2 - ^i) 

Thus the four solutions are 



>i +^2)^ + 4a'] sin'(2^)} 



1/2 



(3.28) 



(3.29) 



u = ±Jp^ + Al{p), 



uj = ±Jp^ + (p). 



(3.30) 
(3.31) 



These again reduce to Mannheim's result if we set /3=0 0. It is interesting to note that these 
solutions are not functions only of 777' + m\ and m\ — m\, as is the case in the relativistic 
regime. 



Fig. ^(b) shows a plot of the dispersion relations for Majorana neutrinos in a medium. 
The dotted and solid curves correspond to the cases with no coupling and with 6=0.2, 
respectively. Again the curves with non-zero coupling are similar to those with no coupling, 
except for the "level repulsion" which occurs in the former case. The parameters in this plot 
are identical to those in the analogous plot for Dirac neutrinos shown in Fig. ^(a). Clearly 
the dispersion relations are quite different in the two cases. 

In all cases examined the minimum of the dispersion relations is always greater than or 
equal to the minimum mass and so again there appears to be no trapping in the Majorana 
case. 

For relativistic neutrinos the exact expressions for the energies may be simplified some- 
what to give 



.:.p-(2/3-a) + |±i 

1 

2p Ap 



{4ap - A2 cos(2^)) + sin2(2^) 
(Aap - cos(20))^ + A^ sm'^{2e) 



nl/2 



1/2 



(3.32) 

(3.33) 



the first of which is in agreement with the analogous expression, Eq. (|3.25| ), for negative 
helicity neutrinos in the Dirac case. 



C. The Vector Model 

In the models considered so far, the fact that the dispersion relations have minima at 
non-zero values of the momentum is due to the chiral nature of the potential in the effective 
Dirac equation. That is, since the potential depends on the spin, the curves are displaced 
to the right or left depending on whether the neutrino's spin is parallel or anti-parallel to 
its momentum. As we have seen, this phenomenon occurs for a single neutrino flavour and 
persists when another flavour is added. It is amusing to note that it is possible to obtain a 
minimum at non-zero momentum even with a purely vector interaction, although this effect 
requires the presence of at least two neutrino fields. Such a model was studied several years 



ago by Chang and Zia ||T2[. The effective Lagrangian is in this case given by 
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-M-aQ7°}V^ = 0, (3.34) 

where the matrix Q is as defined in Eq. ( |3.2| ). The above equation is similar to Eq. ( p.l| ) 
except for the absence of the (1 — 7^) factor which was present in that case. We have also 
set (3=Q for simplicity. The equations governing the dispersion relations may be derived in 
a manner similar to the Dirac case above to yield 



{uj — a cos^ 6^ — — m\ {u — a sin^ 6 



^2 2 
— p — m2 



2o? sin^ Q cos^ Q iu? + — auo + mim2) + sin"^ 9 cos^ 9 (3.35) 



which is independent of the spin and is symmetric under p^—p. A-priori it might then seem 
impossible to generate a minimum at non-zero p. Indeed, for a single neutrino fiavour this 
is the case. For two fiavours, however, something very interesting can happen. Suppose we 
first set 9 to zero and imagine increasing a by so much that the negative energy i/g solution 
overlaps with the positive energy solution. When a non-zero coupling is included, these 
levels repel each other and minima develop near the former crossing points, symmetrically 
placed about the origin. This feature is illustrated in Fig. |^. Thus in this well it 

is possible to have minima in the dispersion relations at non-zero values of the momentum. 
Note that this case is still somewhat different from the chiral cases which we have studied 
above, since the first and second derivatives at the origin are zero and negative, respectively, 
corresponding to a negative effective mass at the origin. This is not the case in chiral 
theories. 



D. Neutrino Oscillations 

We have so far mostly restricted our attention to an investigation of the forms of the 
dispersion relations themselves and have not considered in detail the effects that these would 
have on the oscillations of neutrinos. We first note that in the relativistic regime, the 
standard MSW results are recovered; that is, (i) the neutral current contribution factorizes 
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and (ii) the negative-helicity states obtain the appropriate dispersion relations in matter 
while the positive-helicity states revert to their vacuum dispersion relations. 

For non-relativistic neutrinos, however, the situation is in some sense far more interesting. 
One novel effect which arises in the Dirac case purely as a result of the chiral nature of 
the potential is that in principle one could observe neutrino oscillations with only a single 
neutrino flavour. This could happen since, for non-relativistic neutrinos, the left-handed 
interactions responsible for producing neutrinos would produce both negative- and positive- 
helicity neutrinos. Since these propagate with different phase velocities in the medium, they 
would in general get out of phase with each other, producing oscillations in the probability 
to detect left-handed neutrinos. For relativistic neutrinos this effect disappears since the 
amplitude to produce and detect a positive-helicity neutrino becomes negligible. Note, 
however, that the difference in phase velocities remains and would lead to oscillations if only 
positive-helicity neutrinos could be produced and detected. 

The generahzation of this effect to the two-neutrino case gives the result that in general 
there could be oscillations between four different states. This would lead to an oscillation 
probability which is a superposition of four different oscillation curves. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper we have examined the coherent interactions of a neutrino with a background 
medium by examining the solutions of the effective Dirac equation for the neutrino. A close 
analysis revealed that the dispersion relations corresponding to such a Dirac equation have 
a non-trivial form, even in the simple case in which there is only a single neutrino flavour. 
In particular, we have examined the interesting effects which arise due to the minimum of 
the dispersion relation which occurs for non-zero momentum. We have shown that, quite 
generally for Dirac neutrinos, the minimum value of the energy is less than the neutrino's 
mass, which implies that for any such background there will be trapping of very low energy 
neutrinos. In cases in which the strength of the potential exceeds twice the rest mass of 
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the neutrino in vacuum, neutrino-antineutrino pairs are produced by the electron density 
gradient at the boundary and this affects the reflection amplitude of a low energy neutrino 
in the medium. Our analysis of the case of a single Majorana neutrino flavour showed that 
Major ana neutrinos are not trapped by the medium. 

We have also presented a study of the dispersion relations for Dirac and Majorana neu- 
trinos for the Standard Model with two neutrino flavours. In this case we found that the 
trapping phenomenon persists in the Dirac case but is absent in the Majorana case. The 
neutral current contribution to the oscillation probabilities "factorizes" in the relativistic 
regime, but not in the non-relativistic case. We saw, in fact, that in principle neutrino 
oscillations could occur with only a single flavour of neutrino, due to the different phase 
velocities of the helicity eigenstates. 
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FIG. 1. One-loop diagrams contributing to the neutrino self-energy in the model of Sec. 
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FIG. 2. Dispersion relations for the model considered in Sec. |II A| , with a=3 and m=0, 1, in 
arbitrary units. The dashed and solid curves correspond, respectively, to the m=0 and ttit^O cases. 
Note how the curves in the massive case result from the level repulsion at the level crossings of the 
massless case. 
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FIG. 4. Density-dependent one-loop diagrams contributing to the neutrino self-energy in a 
realistic model with both neutral-current and charged-current couplings. In (a) the "f" stands for 
the contributions due to all fermions in the medium which have neutral current couplings. In (b) 
we have assumed that the medium contains electrons and positrons, but no other charged leptons. 
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FIG. 5. Dispersion relations for two neutrinos in (a) the Dirac case (negative helicity) and (b) 
the Majorana case. In both cases we have set a=1.0, /3=2.5, mi=0.5 and m2=2.0, in arbitrary 
units. The dotted and solid curves correspond to 9=0 and 0.2, respectively. 
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FIG. 6. Dispersion relations in the two-neutrino vector model, with a=3.0, mi=0.5 and 
m2=1.0, in arbitrary units. The dotted and solid curves correspond to ^=0 and 0.2, respectively. 
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